Fast magnetic reconnection in laser-produced plasma bubbles 
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Recent experiments have observed magnetic reconnection in high-energy-density, laser-produced 
plasma bubbles, with reconnection rates observed to be much higher than can be explained by 
classical theory. Based on fully kinetic particle simulations we find that fast reconnection in these 
strongly driven systems can be explained by magnetic flux pile-up at the shoulder of the current 
sheet and subsequent fast reconnection via two- fluid, collisionless mechanisms. In the strong drive 
regime with two- fluid effects, we find that the ultimate reconnection time is insensitive to the nominal 
system Alfven time. 



Magnetic reconnection [l|, Q, the change of mag- 
netic topology in the presence of plasma, is observed 
in space, laboratory, and, most recently, laser-produced 
high-energy-density (HED) plasmas [3|-l6|]. Reconnection 
plays a key role in energy release by plasma instabilities, 
as in solar flares or magnetospheric substorms, and the 
change in topology allows the rapid heat transport associ- 
ated with sawtooth relaxation in magnetic fusion devices 
Q. The observation of reconnection in HED plasmas 
suggests that it may also play a role in inertial confine- 
ment fusion, and indeed recent work has now observed 
the generation of large-scale magnetic fields during iner- 
tial fusion implosions [8]. 

We investigate recent experimental observations 
of fast magnetic reconnection between the HED plasma 
bubbles created by focusing terawatt (TW)-class lasers 
(~kJ/ns) down to sub- millimeter-scale spots on a plastic 
or metal foil. The foil is ionized into hemispherical bub- 
bles that expand supersonically off the surface of the foil. 
Each bubble is found to self-generate a strong magnetic 
field of order megagauss (MG), which forms a toroidal 
ribbon wrapping around the bubble. If multiple bub- 
bles are created at small separation, the bubbles expand 
into one another, and the opposing magnetic fields are 
squeezed together and seen to reconnect (Fig. [1]). The 
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rates of reconnection are observed to be fast^nd unex- 
plained by classical Sweet-Parker theory [3|, Recon- 
nection has also been observed between laser-produced- 
plasma bubbles immersed in a background plasma and 
magnetic field [9], though we focus here on reconnection 
between HED bubbles driven by kJ-class lasers. 

It is of great interest to bring these results in line with 
what is already known about reconnection. As mentioned 
above, there are a number of new features in these laser- 
driven experiments, such as the high energy density in the 
plasma and magnetic field. Perhaps their most notable 
feature is the very strong reconnection drive: the oppos- 
ing magnetic fields are driven together by the expanding 
bubbles at sonic and super- Alfvenic velocities. Such a 
strongly-driven regime has not been previously accessible 
in experiments, but may aid understanding of reconnec- 
tion at the Earth's magnetopause (where the driver of re- 
connection is the super- Alfvenic solar wind) and other as- 
trophysical contexts with high plasma {3 (accretion disks, 
stellar interiors), or with colliding supersonic, magnetized 
flows (supernovae remnants, or at the heliopause). At the 
magnetopause, for example, the strong inflow can drive 
"flux-pileup" reconnection. lol, 



FIG. 1: Magnetic reconnection driven between expanding 
plasma bubbles, view from top. 



Here, we report the results of the first fully kinetic, 
particle-in-cell (PIC) simulations with parameters and 
geometry relevant to these experiments. Our main aim 
has been to understand the apparent Alfvenic reconnec- 
tion rates observed in the experiments. We demonstrate 
that it can be explained as a combination of two effects: 
(1) a strong pile- up of magnetic flux [12] at the shoul- 
der of the current sheet, caused by the sonic inflow and 
corresponding ram pressure that drives substantial com- 
pression of the plasma and magnetic field between the 
two bubbles, and (2) the intervention of collisionless ef- 
fects (Hall current and electron pressure tensor) when 
the current sheet width falls below the ion skin depth 
{di = c/uOpi^ where c is the speed of light, and ujpi the ion 
plasma frequency). Previous discussion of these experi- 
ments has proposed anomalous resistivity [5] , or the Hall 
effect [6], but have not discussed pileup, which we find 
is an essential ingredient. To our knowledge, these laser- 
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driven experiments are the first in which flux pile-up as 
well as collisionless effects play important roles, produc- 
ing highly time-dependent, impulsive reconnection dy- 
namics, and Alfvenic peak reconnection rates. 

In the laser-driven experiments, the reported reconnec- 
tion rates were very fast compared to standard, classical 
predictions from Sweet-Parker theory, and instead were 
consistent with a simple hydrodynamic time based on the 
observed bubble expansion rate. The reconnection in- 
flows (1.3 and 4 x 10^ m/s for Rutherford ^] and Omega 
Q, respectively), ranged from about 0.7 to 3 times the as- 
sociated Alfven speed. (Va = B/^jj^op^ calculated based 
on B fields (100 T, 50 T), and mass density p derived 
from ion species (Al, CH) and electron density (both 
~5 X 10^^ m~^), evaluated at the bubble edge, before 
the onset of the bubble collision.) Despite the difficulty 
making these novel measurements, the overall picture is 
that the reconnection is very strongly driven and occurs 
at Alfvenic to super-Alfvenic rates. Naturally then, the 
inflows are also much larger than allowed by classical 
Sweet-Parker theory, which allows inflows only a frac- 
tion of the Alfven speed, Vsp ~ Va/^S'^/^, where the 
Lundquist number S = LVaPo/t]^ based on system size 
L and classical resistivity 77, roughly 200 and 1000 for 
Rutherford and Omega. Meanwhile, the ion skin depths 
di (40 /im, 40 pm) are equal or larger than Sweet-Parker 
current sheet widths Ssp = L/S^^'^ (15, 30 /im), sug- 
gesting that the experiments are in the regime for two- 
fluid, collisionless reconnection. The two-fluid reconnec- 
tion mechanisms have been extensively observed in simu- 
lation s [13u 14] and have begun to be verified experimen- 
tally [15|. Given this, however, the reported Alfvenic 
reconnection rates remain extraordinary, because even 
these fast, collisionless theories still typically find recon- 
nection inflows "only" near 0.1-0.2 Va- 

We now present 2-d particle- in-cell simulations of bub- 
ble reconnection. The simulations track a pair of ex- 
panding bubbles through their interaction and reconnec- 
tion, which is driven by the substantial energy stored in 
the plasma flow. Our simulations adopt an initial con- 
dition corresponding to a time about halfway through 
the experiments, after the bubbles are created, have ex- 
panded, and have generated their magnetic ribbons, but 



before the pair of bubbles has begun to interact. (It thus 
avoids the physics of the magnetic field generation pro- 
cess, which is a 3-dimensional, Vn x VTg two-fluid effect 
[1, BlO We then model the subsequent expansion and 
interaction of the plasma bubbles. 

We define the following initial conditions. The sys- 
tem here runs, in (x, z) coordinates, over the domain 
[—Lx^Lx] X [—Lz^Lz]. The system is periodic in both x 
and z for computational expediency. It therefore contains 
two half-bubbles, one centered at (0, —Lz) and the other 
at (0, -\-Lz). Define the radius vectors from the center of 
each bubble, r^-*^^ = (x, z + Lz) and r^^^ = {x^z — Lz). 
Then the initial density is + n^^^ + n^^\ where is a 
background density, and the densities n^*^ are 



n^^\x, z) 



'(no-n5)cos2(^) ifr(^)<Ln, 







otherwise. 



(1) 



Here Ln is the initial size of the bubbles. Next, the ve- 
locity field is initialized as the sum of the fields 



_ (Vo sin f if r(^) < Lr 







otherwise. 



(2) 



These simulations adopt uniform initial electron and ion 
temperatures (Tgo and T^o) for simplicity. The magnetic 
field is initialized (Fig. Efa)) as the sum of two toroidal 
ribbons, with 



^^osin(^i%^)r(^)xy 

if r(^) G [Ln-2LB,Lnl (3) 
otherwise. 



Here Bq is the initial strength of the magnetic field, and 
Lb the half- width of the initial ribbons. 

Most parameters can be matched to the experiment, 
except for a heavy electron mass and relatively low speed 
of light, as is typically the case for PIC simulations. 
Lengths are scaled to units of the ion-skin depth to cor- 
rectly retain the two-fluid effects compared to the exper- 
iments. Ln ranges from 20-100 di, for the Rutherford 
and Omega experiments. Lb was found from experimen- 
tal measurements and associated simulations to be 
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FIG. 3: (a,d) Magnetic fields \B\/Bq. (b,e) Reconnection electric fields Ey /VaqBq. (c,f) Generalized Ohm's law Ey — 
-(Vi X B)^ + (l/ne)(j x B)^ - (l/ne)dPe,y / dt - (1/ne) V±Te,^±. Black: sum of Ohm's law RHS; gray: Ey, blue: -Vi,zB:^ 
(with Vi the ion flow); red: {l/ne){j x B)] green: — (1/ne) x dxTe^yx] yellow: — (1/ne) x dzTe^yz] pink: —{l/ne)dPe,y/dt, where 
Pe,y is the momentum density in the electron fluid. The stress tensor T^^ij includes the inertia terms menVe,iVe,j due to electron 
flow along the i and j directions. 



~ 0.1 Ln in Omega. Note that the current sheet width 
is therefore of order the ion skin depth, which is essen- 
tial for the two-fluid reconnection mechanisms discussed 
here. Bq is scaled to the plasma pressure, to achieve the 
high-beta regime of the experiments. 

Particle simulations were conducted with the relativis- 
tic, electromagnetic, explicit particle-in-cell code PScflij. 
We choose the PIC parameters Mi/rrie = 100, T^o = 
Teo = 0.02meC^, and an initial 100 particles per cell. To- 
tal energy is conserved to better than 1% in the simula- 
tions reported. The results here are for collisionless sim- 
ulations. While this is largely consistent with the above 
estimate that reconnection in these experiments is dom- 
inated by collisionless mechanisms, collisions are clearly 
also non-zero, and this will be an important part of future 
modeling as greater experimental fidelity is sought. 

Figure [2] presents the evolution of the magnetic field 
in one simulation with parameters close to those of the 
Rutherford experiments [3], Ln = 20(ii, rib /no = 0.1, 
Ls/Ln = 1/6, Bl/noTeo = 0.25, Vo = 3(TeoM)'/'. 
The image backgrounds show over which are plot- 
ted contours of ^, which is defined such that B(x,z) = 
y x + By{x^ z)y. Panel (a) gives the initial condition 
of Eq. m The next three panels show the evolution at 
later times, with time scaled in units of the dynamical 
time td = Ln/Vo. As the bubbles interact, a reconnec- 
tion X-line forms at the leading point of tangency be- 
tween the bubbles. Near t/td = 0.42, the reconnection 
is about half complete, and the reconnection rate is near 
its maximum. Notice the strong pileup of fiux upstream 
of the reconnection region. In contrast to reconnection 
considered in other contexts, here the strong inflows have 
sufficient force to compress the flux in the current sheet 
by a ratio of about 4. This will be shown momentarily to 
have important consequences for the scaling of the peak 
reconnection rate. Finally, by t/td = 0.76, in (d), re- 



connection is complete, and the two bubbles have almost 
completely merged. Notably, complete reconnection has 
occurred within about 1 dynamical time. 

Figure [3] shows, at two characteristic times (a,b,c: 
t/td = 0.2, (d,e,f): t/td = 0.33), the reconnection 
rate Ey = /dt^ and a detailed cut of how the general- 
ized Ohm's law is fulfilled along a cut across the x-line. 
The time t/td = 0.2 is characteristic of the compression 
or pileup phase: the reconnection rate is small, and the 
electric field is simply supported by the E x B plasma 
flow of the ribbons (c). However, by t/td = 0.33 the re- 
connection is now proceeding at its maximum rate. The 
electric field, of similar magnitude to that earlier in the 
flow, is now sustained (1) in the current sheet by the 
Hall effect (j x B/ne), since the current sheet is now of 
the ion skin depth scale, and (2) at the x-line by the the 
electron stress tensor. Similar to typical results from 2- 
d PIC reconnection studies, within the stress tensor, the 
off-diagonal pressure component {dxTyx)/ne is dominant, 
accounting for 80% of the electric field. 

We find that the reconnection rate /dt sustained 
at this time is super-Alfvenic compared to the nomi- 
nal BqVao^ where B^ is given in the initial condition, 
and Vao is evaluated with no and Bq. Furthermore, 
Ey/ BQVA(d\max — 4.7 is about a factor of 50 higher than 
the typical 0.1 found from reconnection simulations 
including the Hall effect. Here we have normalized to the 
nominal Alfven speeds to correspond to the experiments: 
the magnetic fields measured (and simulated), e.g. in 
Fig. 4(cl) of Ref. [4], are in fact an initial condition be- 
fore the bubbles have begun to interact. 

This point is important because two-fluid reconnection 
mechanisms (the Hall effect and, in these PIC simula- 
tions, the electron pressure tensor), are typically found 
to allow for reconnection inflows of 0.1-0.2 Va- On its 
face, this is not sufficient to explain the apparent Alfvenic 
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FIG. 4: (a) Time for complete reconnection, and (b) peak re- 
connection rate normalized to the instantaneous Alfven speed, 
vs. initial plasma /3e = 2noTeo/5o- 

reconnection observed in the experiments. However, in 
this strongly driven regime, the Alfven speed is also time- 
dependent, due to the compression of the magnetic flux. 
Therefore we also track the instantaneous Alfven speed 
using the maximum upstream magnetic field \B\ and the 
density in the current sheet. Indeed, once one accounts 
for the factor of 4 magnetic flux compression (a factor a 
16 in BVa)^ and the slightly reduced density at the cur- 
rent sheet (n/no — 0.4) one finds that Ey/BVA — 0.25, 
much more in line with typical two-fluid physics results. 
Therefore, we find that previously-established two-fluid 
physics results can also apply to this geometry, once one 
accounts for the flux pile- up. 

Finally, we confirm the importance of flux-pileup in 
this strongly driven regime. Figure |4] shows a sequence 
of simulations in which we have varied the strength of the 
initial magnetic field keeping other parameters fixed. 
The initial inflows become increasingly super-Alfvenic 
compared to the nominal Alfven speed as the initial mag- 
netic field is decreased. We then measure the magnetic 
reconnection rates in two ways. First, in Fig. HJ^a) we 
measure the time for total reconnection, based on the 
time in the simulation to reconnect 95% of the initial 
flux [^0 = (4/7r)L^5o]- Interestingly, we find that the 
total reconnection time only depends quite weakly on Bq ; 
therefore in the strongly driven regime the reconnection 
time is insensitive to the nominal Alfven speed. 

At the same time, we can also study the system at 
the time of peak reconnection. Figure HJ^b) shows the 
peak reconnection rates, normalized to the instantaneous 
Alfven speeds, calculated as above. We find that with 
these corrections, over a wide range of initial magnetic 
fields and therefore flux pileup ratios, the ultimate recon- 
nection rates are consistently about 0.2 BVa-, roughly in 
line with previous two-fluid investigations. 

In past research, flux pileup has been found to op- 
erate when reconnection is driven faster than allowed 
by the reconnection model (resistive MHD [l2|, or 



Hall MHD [Trj ) . If the inflow rates are larger than the re- 
connection rates, then flux must pile up. This is apparent 
in Fig. [31 where the reconnection rate at the x-line is less 
than the flux advection rate from the plasma flow. Mean- 
while, the two-fluid physics result Ey oc BVa oc B^ flnds 
that reconnection rates increase sharply with pileup, so 
that eventually these two effects can come into balance. 
However, pileup can only proceed if no other effects can 
stop the compression process, and an important consid- 
eration here would be the accumulation of plasma and 
magnetic pressure in the current sheet. However, with 
the view of the laser-driven experiments, where the ini- 
tial current sheet densities are much smaller than the 
bubble core densities, we do not believe this would be 
sufficient to prevent reconnection. 

In fact, our simulations find substantial pressure in- 
crease in the current sheet, mostly due to strong heating 
of the current sheet by compression. The observation of 
current sheet heating is qualitatively consistent with the 
Rutherford experimental results, which attributed the 
heating to magnetic energy dissipation (sf. We suggest 
instead that the converging plasma fiows far dominate 
the magnetic fields as an energy source, and that the 
heating is from compression. 

To conclude, recent experiments on magnetic recon- 
nection in laser-produced plasma bubbles are found to 
be in a strongly driven reconnection regime unexplored 
in previous laboratory experiments, where a combination 
of fiux-pileup and two-fiuid effects account for extremely 
fast reconnection. The pile-up is found to be an enormous 
effect, boosting the relevant Alfven speeds to match the 
specified plasma inflow. The strong inflow and flux pileup 
lead to a reconnection time independent of the nominal 
Alfven time. 
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